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ABSTRACT 

We use HST/ACS images and a photometric catalog of the COSMOS field to analyze morphologies 
of the host galaxies of ^400 AGN candidates at redshifts 0.3 < z < 1.0. We compare the AGN hosts 
with a sample of non-active galaxies drawn from the COSMOS field to match the magnitude and 
rcdshift distribution of the AGN hosts. We perform 2-D surface brightness modeling with GALFIT to 
yield host galaxy and nuclear point source magnitudes. X-ray selected AGN host galaxy morphologies 
span a substantial range that peaks between those of early-type, bulge-dominated and late-type, 
disk-dominated systems. We also measure the asymmetry and concentration of the host galaxies. 
Unaccounted for, the nuclear point source can significantly bias results of these measured structural 
parameters, so we subtract the best-fit point source component to obtain images of the underlying host 
galaxies. Our concentration measurements reinforce the findings of our 2-D morphology fits, placing 
X-ray AGN hosts between early- and late-type inactive galaxies. AGN host asymmetry distributions 
are consistent with those of control galaxies. Combined with a lack of excess companion galaxies 
around AGN, the asymmetry distributions indicate that strong interactions are no more prevalent 
among AGN than normal galaxies. In light of recent work, these results suggest that the host galaxies 
of AGN at these X-ray luminosities may be in a transition from disk-dominated to bulge-dominated, 
but that this transition is not typically triggered by major mergers. 

Subject headings: galaxies: active — galaxies: evolution — galaxies: interactions — galaxies: structure 
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1. INTRODUCTION 

Studies over the past decade suggest fundamen- 
tal links between galaxies and their central super- 
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massive black holes (SMBHs). The masses of the 
central SMBHs in nearby galaxies correlate with 
seve ral host bulge properties, including luminos- 
ity (Kormendv fc Richstone 1 119951 : iMcLure fe Dunlopl 
120011: iMarconi fe Hunt I 12003ft. velocity disp ersion 
(iFerrarese fc MerrittJ l2()M . jGebhardt et all 12000 ). and 
mass (Magorr ian et al.l 119981 for a review of these re- 
lations, see Ferrarese 2004). More recent studies extend 
these relationshi ps to galax ies and quasars at redshifts 
from z = 37 (!Ti eti_e^alJ |2004j) to as high as z ~ 4 
(jPeng et alJl2006aU bl). 

These observations indicate co-evolution of SMBH 
mass accretion and host bulge formation processes, 
perhaps through interactions such as AGN feedback. 
Galaxy merger si mulations that includ e a prescription for 
SMBH feedback (|Soringel et al.ll2005h and reproduce the 
MsH-Cbuige relation can recover several observed prop- 
erties of quasar host descendants, including the stellar 
mass function of local elliptical galaxies, and the red 
galax y luminosity function and its variation with red- 
shift (Hopkins et al.l l2006j ). In these models, gas-rich 
mergers drive material toward the central black holes, 
leading to intense star formation and SMBH accretion. 
The nuclear SMBH begins its active phase obscured, but 
feedback energy during the peak accretion phase blows 
away the obscuring material and results in a brief quasar 
phase. The blowout phase coincides with a rapid trun- 
cation of star formation throughout the host galaxy. Re- 
cent observational studies lend credibility to this picture 
by hinting tha t star-formation quenching coincides with 
AGN activity (iBundv e t alJl2007t [ Silverman et"all l2008t 
iTremonti et al.ll2007ft . ~ 

This picture of SMBH-host co-evolution relies on a hy- 
pothesized merger mechanism for fueling active black 
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holes, and incorporates predictions that gas-rich disk 
galaxies merge to form luminous starbursts, eventually 
evolving into massive elliptical galaxies. Other possi- 
ble fueling mechanisms includ e less-direct tidal inter- 
actions with nearby galaxies (|Menci et all I2004T ) and 
instabilities in a quiescent galaxy's gaseous disk (e.g. 
iHopkins fc Hernquist I [20061 . In this study, we explore 
the interaction mechanism for intermediate-luminosity 
AGN in the COSMOS field by analyzing their environ- 
ments and host galaxies. 

Previous similar work e mploying HST survey data has 
found conflicting evidence. iGrogin et alJ ([2005T I analyzed 
^100 X-ray selected AGN at redshifts up to ~1.3 in 
the GOODS fields and found no significant differences 
between their host structural properties and companion 
fractions and those measured for matched control sample 
galaxies. iPierce et all (|2007l ). on the other hand, exam- 
ined ~60 X-ray and IR-selected AGN with 0.2 < z < 1.2 
in the Extended Groth Strip and found that AGN are 
marginally more likely than control galaxies to have 
nearby companions. In both studies X-ray selected AGN 
are found to reside predominantly in host galaxies with 
bulge-dominated morphologies, generally in agreement 
with work at low and hi gh redshift qu asar host stud- 
ies (jJahnke et al l l2004bl : ISanchez et all I2004D . Inves- 
tigating larger-scale environ ments of 5 2 quas ars from 
the DEEP2 redshift survey, ICoil et all (|2007h showed 
that the quasar-galaxy cross-correlation function at z ~1 
closely resembles the galaxy auto-correlation function at 
all scales, and that the relative quasar bias traces that 
of blue galaxies better than red galaxies. This might 
suggest that high luminosity AGN reside in blue bulges 
(jJahnke et alJ I2004U ISilverman et alJ 12008ft that have 
not yet migrated to the high-density environments typi- 
cally found for massive red galaxies. 

In the nearby universe, where precise techniques al- 
low detailed studies of tens of thousands of AGN hosts 
(with, e.g., the Sloan Digital Sky Survey), active nu- 
clei reside in the most massive galaxies, with structural 
properties similar to early- type galaxies but with rela - 
tively young stellar populations (|Kauffmann et aLll2003f ). 
Examining imag es of ~100 of the most luminous AGN 
within z < 0.1, Kauffman n et al.l found that they oc- 
cupy roughly equal fractions of blue spheroids, single disk 
galaxies, and disturbed/interacting galaxies. Further- 
more, low-redshift close galaxy pairs that exhibit strong 
indications of interaction are more likely to include AGN 
than pairs without interaction indicators (jAlonso et al.l 
2007j). At the same time, black hole accretion activity is 
significant ly larger for AGN w ith bright companions than 
otherwise (jAlonso et al.ll2007h . and quasars are found to 
have local galaxy overdensities withi n 100 kpc in excess 
of that seen in non-active g alaxies (ISerber et alJ [20061 ) 
and lower-luminosity AGN ([Strand et alJl2008f ). Com- 
bined with earlier imaging studies showing that a mi- 
nori ty of local Type 1 AGN are undergoing interactions 
(e.g. LDe Robertis et alJfl998f >. these studies suggest that 
mergers and interactions play some role in fueling AGN, 
but not necessarily a dominant one. 

These low-redshift AGN may represent a different pop- 
ulation than that found closer to the peak of AGN ac- 
tivity, z ~ 2, with different fueling mechanisms in play. 
The typical AGN in a local survey like the SDSS is less 
luminous and possibly hosted by a galaxy in a differ- 



ent evolutionary state than AGN selected at higher red- 
shift. Intermediate- and high-redshift AGN (z > 0.5) 
are typically found in bluer, more extended, and more ir- 
regular galaxies than their low-redshift c ounterparts (cf. 
Pahnke et al"1l2004bl : [Sanchez et alJl2004f >. By using sam- 
ples with a wider redshift range we can constrain the 
dominant mechanism in the history of AGN fueling, and 
potentially uncover its evolution over cosmic time. Sam- 
ples at z > 0.3, however, remain limited to a few dozen 
objects selected through a variety of methods (see stud- 
ies mentioned above). In the present study, we extend 
this moderate-redshift sample and compare the effects of 
different selection techniques. 

Using the extensive data of the Cosmic Evolution 
Survey (COSMOS), we determine basic properties of a 
sample of AGN hosts to probe their co-cvolution with 
SMBHs. Selection of the AGN sample and data used 
for the analysis is described in fj2j The analyses of the 
host galaxies are presented in §Sj3]and[H followed by a 
brief discussion and conclusions. For cosmological cal- 
culations, we adopt h = 0.75, (where Hq — 100ft. km 
s _1 Mpc -1 is the Hubble parameter), £Im = 0.3 (matter 
density parameter), and Q\ = 0.7 (cosmological constant 
density parameter). 

2. DATA AND SAMPLE 

COSMOS (jScoville et aUl2007h . a Hubble Space Tele- 
scope Treasury project, includes coverage of a 2 square 
degree field from X-ray wavelengths to UV, optical, IR, 
and radio. The cornerstone data set, which we use for 
the bulk of our analysis, consists of 583 orbits taken 
with Hubble's Advan ced Camera for Surveys (ACS) with 
the F814W filter (see iKoekemoer et al.ll2007l for a com- 
plete description). Ancil lary observations incl ude XMM- 
Newton X- ray imaging (Hasi nger et aT][200l . VLA ra- 
dio maps (ISchinnerer et all l2007f)~and VLT/VIMOS 
([Lilly et alJ 120071) and Mage llan/IMACS optical spec- 
troscopy ([Trump et all I2007L Trump et al. 2008, in 
preparation) . 

Our sample selection focuses on AGN candidates in 
the COSMOS field with spectroscopic redshifts. An ob- 
ject is identified as an AGN candidate through detec- 
tion as an X-ray point source above the ~ 10 -15 erg 
cm~ 2 s " 1 flux limit in the 0- 5 -2 keV or 2-10 ke V flux 
bands ([Cappelluti et all 120071 iBrusa et all I2007H . or a 
radio source aboye the 0.1 mJy flux limit at 1.4 GHz 
(Schi nnerer et al.l 12007). Optical counterparts to these 
candidates with Iab < 24 are followed up in the Mag- 
ellan/IMACS spectroscop ic survey, whose firs t season of 
observations is detailed in lTrump et a l. (2007), including 
emission line identification and redshift determination for 
^350 objects. We include additional sources from the 
second season of IMACS observations, as well as compan- 
ion observations using MMT/Hectospec (Trump et al., 
in preparation). Targets with successful redshift deter- 
minations are separated into 4 primary spectral classes: 
1) broad emission line AGN; 2) narrow emission line ob- 
jects; 3) red galaxies, with detectable continua but no 
emission lines; and 4) hybrids showing narrow emission 
lines superposed on a red galaxy continuum. 

From the spectroscopically-confirmed AGN candi- 
dates, we restrict our redshift range of interest to 0.3 < 
z < 1.0. The upper limit arises from our limited abil- 
ity to adequately analyze targets at high redshift us- 



AGN Morphologies 3 



TABLE 1 
Numbers of AGN in the samples. 



AGN sample No. 


of objects a 


w/ ACS images 


Successful 2-D fits 


AGN candidates w/ spectra 


~1300 


*** 


*** 


Spec, redshifts 0.3 < z < 1.0 


459 


391 


314 


Type 1 


48 


36 


19 


X-ray Class 2 


83 (95) 


73 


57 


X-ray Class 3 


48 (62) 


11 


36 


Radio Class 2 


40 (120) 


32 


30 


Radio Class 3 


82 (134) 


71 


62 



a Numbers in parentheses include those objects in the X-ray sample with log(Lx) < 42, 
and those objects in the radio sample not classified as AGN. Objects failing to satisfy 
these AGN classifications are excluded from the other columns, as well as the analyses 
in this work. 



ing single-orbit ACS data, and because the ACS F814W 
bandpass shifts into the rest-frame UV at z > 1, bias- 
ing morphological characterization. Based on our sim- 
ulations described in section 3.1, typical AGN hosts at 
this redshift have recovered magnitude uncertainties of 
~0.4 magnitudes. The lower redshift limit is a practi- 
cal one, applied to limit ourselves to objects which have 
a large number of corresponding inactive galaxies and 
whose environments can be analyzed adequately within 
a 2 square degree field. Additionally, some of the ancil- 
lary survey boundaries extend beyond those of the ACS 
observations, so ^60 AGN candidates were identified for 
which we have no Hubble data for host analysis. These 
objects are excluded from the sample. 

Determining which candidates truly host an AGN is 
non-trivial. Those objects with broad emission lines are 
the easiest to classify as AGN, but we treat them sepa- 
rately due to the uncertainties in analysis of their hosts 
(see below). A significant fraction of the narrow- line ob- 
jects may be star-forming galaxies rather than genuine 
AGN (Trump, private communication), and the IMACS 
spectra lack sufficient signal-to-noise and spectral range 
to make suc h a distinction using line-ratio diagnostic di- 
agrams (cf. iBaldwin et ailll981l:IVeilleux fc Osterbrockl 
119871 iKewlev et al.ll200H iKauffmann et alj|2003f) . Fur- 
thcrmore, many X-r ay selected candid ates, which are ex- 
pected to be AGN (jMushotzkv II2004D . exhibit no emis- 
sion lines, falling into the red galaxy spectral classifica- 
tion. This might occur due to obscuration of the regions 
emitting optical spectral lines, misidentification of the 
optical counterparts, or misplacement of the slit when 
performing optical spectroscopy. 

With these complications in mind, we present our find- 
ings separately for various sub-samples, and we use the 
following terms to describe them. "AGN candidates" 
refers to the full sample of candidate AGN with spec- 
troscopic redshifts, and includes 459 objects. "Class 1" 
or Type 1 "broad line" AGN refers to those candidates 
with broad optical emission lines, easily distinguishable 
as AGN or quasars, and including 48 objects. "Class 
2" or "narrow emission line" objects refers to AGN can- 
didates with narrow emission lines in their spectra, in- 
cluding 120 radio-selected objects and 95 X-ray-selected 
objects. This category generally includes both Type 2 
AGN as well as star-forming galaxies. "Class 3" or "non- 
emission line" candidates are those whose spectra look 
like red galaxies with no emission lines, including 134 
radio-selected and 62 X-ray selected objects. In cases 



where an object has been detected both in radio and X- 
ray emission, we include it in the X-ray category. These 
sub-samples are summarized in Table [TJ 

As in pre vious AGN stud i es with similar red- 
shift ranges (iPierce et alJ l2007t ISilverman et all 120081 : 
iBundv et al.l 120071 ). we adopt a cutoff in X-ray luminos- 
ity above which ob jects are likely to be AGN-dominated 
(jBauer et al.ll2004D . After excluding broad-line AGN, we 
take those candidates with L 2 -iokeV > 10 42 er g s s_1 as 
the most likely to harbor accreting black holes. Here, 
^2-iofceV is the X-ray luminosity in the 2-10 keV band. 
We estimate this luminosity for all X-ray sources by con- 
verting the observed XMM-Newton X-ray fluxes via a 
k-correction. To do so, we assume an X-ray power law 
slope of r = 1.9 and perform the k-corrections using 
spectroscopic redshifts. Of the 95 X-ray-selected Class 2 
objects, 83 satisfy the luminosity cut given above, and 
we refer to them as "X-ray Class 2" AGN. Of the 62 
X-ray-selected class 3 objects, 48 satisfy the luminosity 
cut, and we refer to them as "X-ray Class 3" objects. In 
some parts of our analysis, we combine these two sam- 
ples together, and generally refer to them as X-ray AGN. 
Those objects which do not satisfy the luminosity cut 
are excluded from the following analyses. We note here 
that some X-ray point sources have more than one possi- 
ble optical counterpart, with the most likely counterpart 
chosen. This ambiguity applies only to three objects in 
our final sample, so we suspect it has little effect. 

For the radio-selected AGN candidates, we use the 
novel technique for s eparating AGN and star-forming 
galaxies described in (|Smolcic et al.ll2b08f) . Briefly, the 
technique uses a combination of morphology and rest- 
frame colors of optical counterparts to the radio sources 
to classify them as QSOs, stars, star-forming galax- 
ies, AGN or high-redshift galaxies. From the mod i- 
fied Stromgre n photo metric system (|Odell et al.l [2002) , 
ISmolcic et alJ ()2006l ) derive two principle-component 
combinations of rest-frame colors which "optimally quan- 
tify the distri bution of galaxies in the rest-frame color- 
color space" (|Smolcic et alJ I2008I ). Because the emis- 
sion line strengths correlate with galaxy spectral en- 
ergy distributions, these color parameters can mimic line 
diagnostic diagrams' ability to separate emission line 
galaxies into AGN, c omposite, and star-forming objects. 
ISmolcic et alJ (|2008l) calibrate this color classification 
scheme using galaxies from SP SS matched to the N RAO 
VLA Sky Survey at 1.4 GHz (ICondon et alJll998h. and 
apply it to the COSMOS radio sources (|Schinnerer et al.l 
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2007). The low-redshift calibration suggests that the 
sample classified as AGN contatins ~5% star-forming 
galaxies, ~15% composite objects, and ~80% AGN, and 
these AGN comprise ~90% of the total population of ra- 
dio AGN. When we apply the classification scheme to 
the radio-selected objects in our sample, we find that 
122 are identified as AGN, including 40 Class 2 objects 
and 82 Class 3 objects. We refer to these as radio AGN 
throughout the paper. Roughly 30% of our radio-selected 
AGN candidates are not classified eith er as AGN or star- 
forming because ISmolcic et a l. ( 2008) a) use a conserva- 
tive search radius of 0.5 arcseconds when identifying opti- 
cal counterparts, whereas the spectroscopic follow-up in- 
cluded objects as far as 1 arcsecond from the peak radio 
emission, and b) they exclude from their radio-optical 
sample a fraction of objects (15%) that are photometri- 
cally flBgg_edJn_the_COSMOS photometric redshift cat- 
alog (CapakeFlil] {2007, ~ 30% of which have availabl e 
spectroscopic redshifts; see Tab. 1 in lSmolcic etHl2 008) . 
The remaining objects, mostly Class 2, are classified as 
star-forming galaxies. 

We include the broad line AGN for comparison only. 
Because the redshift distribution of broad line AGN 
peaks at higher redshift than the Type 2 distribution, 
the Type 1 sample includes only ~50 objects. Further- 
more, the nuclear component in broad-line AGN tends 
to dominate the total optical flux more than in narrow- 
line AGN, so some of the host-fitting techniques that we 
apply suffer from more serious systematic uncertainties. 
These factors hamper our ability to compare with confi- 
dence the structural parameters and interaction indica- 
tors between Type 1 and Type 2 AGN hosts. We apply 
the present analyses to all objects for completeness. 
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Fig. 1. — Top panel: redshift distributions of narrow-line (solid 
line) and broad-line (dotted line) AGN cadidates, along with can- 
didates with no emission lines (dashed line), from seasons 1 and 
2 of IMACS observations. Bottom panel: redshift distributions of 
narrow line AGN candidates selected using radio emission (dashed 
line) and X-ray emission (solid line). 

Our AGN sample thus includes ^200 objects with 
spectroscopically confirmed redshifts in the range 0.3 < 
z < 1.0, narrow emission line identifications, and ACS 
imaging. Figure [T] shows the redshift distributions of 
AGN in our samples, extending to lower redshift for ref- 
erence. We suspect the redshift peaks near z = 0.3 
and z = 0.7 are assoc iated with large sca le structures 
in the COSMOS field (|Scoville et al.ll2007t) . Broad-line 
objects are plotted for comparison, and the distributions 



of radio-selected and X-ray-selected narrow-line objects 
are compared. 
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Fig. 2.— Example Hubble/ACS F814W images of AGN with 
0.3 < z < 1.0. The bar in the lower left corner of each panel is 
one arcsecond long (~ 5h~ physical kpc), and the pixel brightness 
scale is logarithmic. We show the redshift (z) in the lower right 
corner. 

We apply several techniques, described in the fol- 
lowing sections, to the ACS images to analyze the 
AGN host galaxies. Each image combines data from 
4 sub-exposures, achieving a 5 a point source sensitiv- 
ity of ~27.2 magnitude s in Iab with the F814W filter 
(|Koekemoer et all [2007). We use the original drizzled 
images with 0.03 arcsecond per pixel scale, not rotated to 
north-up orientation, with multidrizzle parameters cho- 
sen to best preserve the original point spread function 
(PSF). Figure [2] shows example cutout images of AGN 
in our sample. To complement the morphological anal- 
ysis, we use the COSMOS photometric re dshift catalog 
(jMobasher et all 120071 : ICapak et ai1l2007D as a filter to 
search for near neighbors to the AGN. 

In order to compare the AGN sample with non-active 
galaxies, we also identify a control sample using the COS- 
MOS photometric redshift catalog. Within a survey as 
large and deep as COSMOS, an ideal control sample 
would include galaxies with the same mass and redshift 
distributions as the AGN. Due to possible contamination 
from the nuclear point sources, we cannot use the pho- 
tometric information to determine mass-to-light ratios, 
and thus masses, for our AGN candidates. Therefore, we 
use the best-fit apparent magnitude of the AGN hosts 
to match the luminosity and redshift distribution of the 
control sample so that k-corrections can be ignored. For 
each AGN in the sample, we find ten non-active galaxies 
with similar redshifts (Az < a z , where a z is the error in 
the photometric redshift for the control galaxies) and ap- 
parent magnitudes (AI AB < 0.3, using the ACS F81AW 
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detections) to those of the AGN host. We perform anal- 
yses on the control galaxies in a comparable fashion to 
the hosts, as described in the following sections. For 
some comparisons, we find it illustrative to separate the 
control sample into early and late spectral types using 
the photometric redshift catalog T p hot parameter, which 
classifies galaxies on a scale from 1.0 (red elliptical) to 
6.0 (starburst) using photometric fits to the galaxy spec- 
tral energy distributions. We divide the control sample 
at Tphot = 2.0, corresponding to an Sa/Sb spectral type. 

3. MORPHOLOGICAL ANALYSIS 

Using the deep, high-resolution ACS COSMOS images 
we attempt to determine properties of the host galax- 
ies of our AGN sample. The high angular resolution of 
Hubble's diffraction-limited imaging allows us to sepa- 
rate host galaxy light from that of the nucleus. Using 
only the I-band images, we can constrain the magni- 
tude, scale, radial light profile, and orientation of the 
host galaxy. 

In this work we use 2-d imensional surfac e bright- 
ness fitting (with GALFIT, iPeng et all l2002f) to mea- 
sure AGN host properties. To understand systematic 
uncertainties in the surface brightness fitting, we simu- 
late AGN images and apply identical fitting techniques. 
After decomposing the images into AGN point source 
and galaxy light, we measure the asymmetry and con- 
centration of the underlying host galaxies and compare 
them directly with the non-active control galaxies. We 
first describe our simulated AGN images, then explain 
the techniques we use for 2-D surface brightness fitting, 
and describe the results of our 2-D fits. Finally, we dis- 
cuss the asymmetry and concentration measurements. 

3.1. Simulations 

We performed two types of simulations to help under- 
stand systematic uncertainties in our analysis. One suite 
of simulations aims to quantify our ability to reliably re- 
cover parameters in 2-D surface brightness models. The 
other suite examines the effect of point-spread function 
(PSF) variation and mis-application on the analysis. In 
this way we isolate the impacts of the two most impor- 
tant problems with 2-D fitting. 

3.1.1. PSF Variations 

In performing 2-D fits to galaxy images, a PSF must 
be supplied to convolve with the galaxy model image. 
Fits of AGN images are especially sensitive to the PSF 
due to the sometimes bright, nuclear point source whose 
light is superimposed on that of the host galaxy. The 
ACS instrument's PSF ellipticity and size are known to 
vary both tempo rally and across th e CCD at the level 
of a few percent (jRhodes et al. 2007). Our solution, de- 
scribed in 33.21 includes sets of model PSF grids. Be- 
cause systematic uncertainties in the PSF can dominate 
the morphological classification of compact AGN hosts, 
we have performed a series of simulations to test how our 
ability to recover host properties varies with PSF. 

We simulate AGN images by superimposing a real star 
extracted from an ACS image onto sets of simple model 
galaxy images with varying parameters whose ranges are 
similar to those of the AGN sample. The model galaxy 
images are created using GALFIT, with effective radius 



and magnitude randomly chosen from uniform distribu- 
tions in the ranges 0.15" < r e ff < 2.5" and 19 < m < 
24, respectively. GALFIT convolves the specified galaxy 
model with chosen stellar PSF to yield the model galaxy 
image. The star image is then scaled to a random mag- 
nitude with 16 < m < 25 and added to the galaxy model 
image. We create 2000 such simulated images with expo- 
nential disk profiles, and an additional 2000 with deVau- 
couleurs profiles. We use four different real star images 
for the PSFs, with 500 simulated images per star. One 
of the four stars was chosen to be near the limit of satu- 
ration, and results for the corresponding 500 simulations 
are obviously skewed and thus ignored in later discussion. 
We refer to the simulations as "PSF simulations," and 
we use them below to characterize the systematic effects 
of PSF variation on the results of 2-D surface brightness 
fits. To check whether the simulations created with these 
four stars adequately encompass the full range of PSF 
variations, we also created 500 simulated AGN images 
using 50 different stars (10 simulated images per star), 
each taken from a different ACS tile and a different de- 
tector position. As the results of fits to these images 
closely mimic the results obtained with the original four 
stars, we leave them out of the discussion below. 

3.1.2. Parameter Recovery 

Even if we have applied a perfect PSF, signal-to-noise 
limits our ability to recover fit parameters accurately, and 
parameter uncertainties are dominated by systematic ef- 
fects. To gauge the robustness of recovered parameters 
and assign appropriate uncertainties, we have performed 
a set of simulations where the PSF remains constant. 

We simulate 2000 AGN images with exponential disk 
hosts, and another 2000 with deVaucouleurs hosts, with 
the same range of parameters as for the PSF simulations. 
To better represent the background sky noise, which is 
the dominant noise component in our images, we ran- 
domly add cutouts from COSMOS images. While these 
background images will sometimes include contaminat- 
ing galaxies, the same is true for our real AGN images 
and the overall effect of the galaxies is minimal. We refer 
to these simulations as "Recovery Simulations," and by 
performing 2-D fits on them we characterize the uncer- 
tainties in our best-fit AGN parameters due to noise. 

3.2. 2-D Surface Brightness Fitting 

We use GALFIT (jPeng et al.1 l2002f l to fit models to 
AGN images in the sample. For each image, we model 
the nuclear point source as a point spread fu nction, and 
the h ost galaxy as a single Sersic function (see lPeng et al.l 
2002, for details of the functional forms of different mod- 
els in GALFIT; Sersic 1968). In short, the Sersic func- 
tion is a general galaxy model which encompasses a range 
of more specific models through the variation of an in- 
dex, n. The Sersic function with n = 1 is equivalent 
to an exponential disk model, whereas a Sersic func- 
tion with n = 4 is equivalent to the de Vaucouleurs 
(r 1 / 4 ; Ide Vaucouleurs fc Capaccioli Ill979f ) profile, which 
describes typical galactic bulges and early-type galaxies. 
The fit results include point source position and magni- 
tude (m p ), along with the host galaxy magnitude (ro^), 
effective radius {rh), Sersic index (n), axis ratio (b/a), 
and position angle in the image. Because some of our 
AGN candidates may not have a nuclear point source, 
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we also performed fits which excluded the point source 
component and used just a single Sersic galaxy model. 

Running GALFIT requires an initial guess of each of 
the best-fit parameters, an input image, a point-spread 
function image, and a sigma image. Input AGN images 
are cut directly from the original ACS images, with a 
cutout image size corresponding to 35 hr 1 kpc comov- 
ing (~17" for z = 0.3 and 6" for z = 1.0; larger and 
smaller image sizes were attempted, with no impact on 
the resulting best-fit parameters). In order to generate 
initial guesses in an automated way, we developed a pro- 
ce dure similar to that u sed in GALAPAGOS, described 
bvlHaufiler et all (120071) First we run Source Extractor 
(|Bertin fc Arnouts 1119961 ) on the cutout image. For every 
extracted source, we generate an elliptical mask image 
using the Source Extractor FLUX_RADIUS, ELONGATION, 
and THETA_IMAGE output parameters for that source. To 
ensure conservative estimates of galaxy boundaries, we 
set the mask semi-major axis to 2xFLUX_RADIUS. Since 
the cutouts are centered on the AGN coordinates (which 
are taken as the optical counterparts of X-ray sources), 
we select the extracted source nearest the center of the 
image as the target AGN. For each additional source in 
the image, we include it in the 2-D fit if and only if its 
mask overlaps the mask of the AGN, and otherwise we 
simply mask it out of the image. All added objects are 
modeled as single Sersic function profiles. Finally, we 
identify the brightest pixel within the AGN mask as an 
initial guess for the location of the point-spread function 
component. 

For Sersic function profiles included in a fit, we 
estimate the initial parameters using results from 
Source Extractor. The effective radius is set to 
r h = 0. 162xFLUX_RAD i US ba sed on the simulation re- 
sults of Haufi ier et al.l (|2007f ). Magnitude guesses are 
set to MAG_BEST, the axis ratio determined from the 
ELONGATION parameter, and the position angle computed 
from THETA_IMAGE. We constrain the Sersic index to lie 
between 0.5 and 8, the magnitude to stray no further 
than 5 magnitudes from the initial guess, and the effec- 
tive radius to be less than 500 pixels (15 arcseconds). For 
the AGN host Sersic component, we constrain the effec- 
tive radius to be less than half the image width. For the 
PSF component, we set the initial magnitude to 3 mag- 
nitudes fainter than the AGN host component (based on 
typical previous fits of the AGN sample). We constrain 
the PSF magnitude to within 10 magnitudes of its initial 
value, and the position to lie within 5 pixels of its initial 
location. The sky value for the image is held constant 
based on the sky subtraction of the original COSMOS 
ACS images. We tested several methods for comput- 
ing the sigma image, including conversion of the weight 
images output by MultiDrizzlc which correspond to the 
ACS tiles, as well as an empirical determination of the 
noise based on the rms signal of regions of sky around 
each AGN candidate. Differences in the choice of sigma 
images lead to uncertainties which are small compared 
to those introduced by PSF mismatch and other effects 
described below. 

We choose the parameter constraints largely by con- 
vention, but also to ensure that they fully encompass the 
reasonable ranges of the parameters. As described be- 
low, we exclude from further analysis those fits which run 
into constraints, since these did not find a true best-fit 



and the parameters are likely unphysical. We attempted 
some variations on the fitting constraints; notably, we 
performed fits without constraints and fits where we re- 
stricted the Sersic index to n < 5 rather than 8. Fits 
without constraints fail to converge with a higher fre- 
quency than those with constraints, although this oc- 
curs mostly because our model poorly matches the real 
light distribution in cases of failure. Fits with a more 
restricted Sersic index yield comparable results to those 
obtained with the original n < 8 constraint. A vast ma- 
jority of fits with n > 5 in the original fits (~80 objects, 
including 25 which run into the n = 8 constraint) run into 
the constraint when we restrict n < 5. Placing the con- 
straint at n < 5 would effectively eliminate those objects 
from our further analysis. However, other parameters 
of the fits (e.g. host magnitude) may yield reasonable 
and useful estimates even when a fit runs into the con- 
straints, and these other parameters can be sensitive to 
the constraints chosen. We find that the more restric- 
tive Sersic index constraint yields host magnitude esti- 
mates systematically 0.13 magnitudes higher (dimmer) 
than the original constraint, with a scatter ~0.5 mag- 
ni tudes. This syste matic offset differs from the results 
of iKim et~ al. (2008) because the objects here are heavily 
skewed toward host galaxy-dominated images. Further- 
more, radius estimates in the restrictive Sersic index case 
are a median of 15 pixels (0.45 arcseconds) smaller than 
in the original case. These relatively minor differences 
do not affect our main conclusions. 




Radius (pixels) 



Fig. 3. — Normalized radial point spread function profiles. The 
shaded region (upper panel) and dashed lines (both panels) show 
the variation from the 10th percentile to the 90th percentile of 
profiles for 50 stars extracted from ACS images. Thick solid lines 
show the profiles for three TinyTim PSF models. In the lower 
panel, we divide each profile by the median profile (dotted straight 
line at 1.0) of the 50 real stars. 

For the point-spread function (PSF) solution to the 
ACS imaging , we adopt the PSF grids described in 
iRhodes et alT(|2007f l. These authors were motivated by 
the demands of detecting weak lensing signals, which re- 
quire characterization of image cllipticiti es at the ~ 1% 
level. Briefly, they use TinyTim software (jKrist Il2003f ) to 
create a PSF model at each of ^4000 points in a regular 
grid, and develop several such grids corresponding to dif- 
ferent focus offsets of the Hubble Space Telescope during 
exposure. For each COSMOS ACS image, the best-fit fo- 
cus position is obtained by simultaneously matching the 
shapes of model PSFs to ~10 bright stars chosen from 
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the image. This process is found to be repeatable to 
an accuracy of ~1 /jm in focus position. Figure [3] ex- 
hibits the variation in the PSF profiles for 50 real stars 
selected from different ACS images at different detector 
positions, along with three TinyTim PSF profiles. The 
TinyTim PSFs can generally encompass the variations 
of the real PSFs, although at large radii they system- 
atically underestimate the flux level of real PSFs. For 
an AGN at any position in an image we use the nearest 
model PSF from the best-fit grid. The results of 2-D fit- 
ting of simulations described above exhibit some of the 
resulting systematic effects of inappropriate PSF choice. 
These effects are described below in §3.31 

After initially fitting all the AGN candidates with 
GALF1T, we determined whether the results had run 
into the boundaries set by the parameter constraints. 
Many of those objects for which this is the case have 
compact light profiles, so we adjusted the initial guess 
file such that the point source magnitude equals the host 
galaxy magnitude and the host radius is double the orig- 
inal guess, and ran GALF1T once again on those objects. 
Finally, we visually inspected all of the resulting model 
images and residual images, subjectively assessing the 
quality of the fit, and in some cases attempting to rem- 
edy a failed fit. This typically entailed masking out a 
nearby star or galaxy whose light was contaminating the 
image beyond its original mask. We discuss fit results in 
the following section. 

In order to place constraints on those objects for which 
the fits failed altogether, we used a simple point source 
subtraction method. First, we fit each AGN with a sin- 
gle point source component and no galaxy component in 
GALF1T. Then we subtracted the best-fit point source 
from the image. On the residual image, we identified 
pixels whose flux values changed from positive on the 
original image to negative after subtraction (indicating 
over-subtraction), and set those pixels to zero flux. Then 
we estimated a lower limit for the host galaxy magnitude 
by using aperture photometry to measure the residual 
flux in an aperture with a two arcsecond diameter. 

We followed similar procedures with the sample of con- 
trol galaxies as with the AGN candidates themselves. We 
first fit the galaxies without a central point source com- 
ponent. To mimic the process of fitting AGN, we then 
superimposed a point source and applied a fit procedure 
identical to the one used for the AGN candidates. Since 
each AGN candidate has ten matched control galaxies, 
each control galaxy is matched to a particular AGN can- 
didate. We thus determined the brightness of the super- 
imposed point source by using the best-fit point source 
magnitude from the AGN candidate fit. Thus, the fitting 
performed on the control galaxies is well-matched to that 
performed on the AGN candidates. 

3.3. 2-D Fitting Results 

Since the formal statistical uncertainties output by 
GALF1T tend to underestimat e the true uncertainties, 
we follow HauBlcr et alj (|2007l ) and use the mean sur- 
face brightness as a proxy for image signal-to-noise to 
diagnose the reliability of recovered fit parameters. The 
mean surface brightness is defined here as /i = rrih + 
2.5\og(2irrf l b/a). Figure 0] shows the distributions of 
measured mean surface brightness for both our sample 
of AGN and for matched control galaxies. The distribu- 



tions have a mean near 22 magnitudes per square arc- 
second, with a standard deviation of approximately 2 
magnitudes per square arcsecond. To connect to more 
physically meaningful galaxy characteristics, we use mea- 
surements of our large sample of control galaxies. A typ- 
ical ~ L* elliptical galaxy in the middle of our redshift 
range at z = 0.7 has an effective radius ~ 0.5 arcsec- 
onds (~ 2.5/i~ 1 kpc) and \x ~ 20.5 magnitudes per arc- 
second, while an exponential disk galaxy has r/, ~ 0.75 
arcseconds (~ 3.7/i _1 kpc) and // ~ 22.1 magnitudes 
per arcsecond. In general, however, the reliability of a 
given best-fit parameter measurement depends in a com- 
plicated way on the other parameters in the fit. 

In particular, the inclusion of a central point source in 
the fit model increases the complexity of the relationships 
between best-fit parameters, so we consider the bright- 
ness of the galaxy component relative to the point source 
component as another important diagnostic. We show 
this as a difference in magnitudes in the right panel of 
Figure 21 plotting the distribution of m p — rrih (point 
source magnitude minus host galaxy magnitude for the 
AGN candidates in the sample. High values of m p — rrih 
correspond to host-dominated images, and low values 
correspond to point-source dominated images. A large 
majority of the AGN in our sample are dominated by 
their host galaxies. This begs the question: should we 
really include the point source component of the fit at 
all, or instead use just a single galaxy component? We 
address the issue by fitting a sub-sample of ~500 con- 
trol galaxies with a galaxy + PSF model. By including 
a point source component in the fit, we can determine 
the relative point source flux at which our fitting pro- 
cedure spuriously identifies a PSF. The distribution of 
spuriously recovered m p — rrih peaks near 4.5, with a 
broad range from ~2 to 6 (dashed line in right panel of 
Figure 0]) ■ We show the cumulative distribution for all 
control galaxies, early-type, and late-type galaxies in Fig- 
ure About 16% (corresponding to the 1-sigma bound- 
ary) of these fits have m p — rrih < 3.4, and ~7% have 
nip — rrih < 3.0. We establish these levels as limits to 
our ability to recover real point sources in AGN images. 
If any given AGN fit has m p — rrih > 3.4 (or 3.0, more 
conservatively), then it is consistent with normal galax- 
ies, lacking a real point source. We take those AGN fits 
with m p — rrih < 3.4 (3.0) to have a real point source 
detection. 

Formally, we can statistically determine the probability 
that an AGN fit with m p — rrih < L includes a real point 
source detection. If we let R represent a real point source 
detection, and P represent a positive measurement that 
m p — rrih < L, then we want to determine p(R\P), the 
probability that an object has a real point source given 
a positive measurement. From Bayes' Theorem, this is 

p{R]p) _ EfflW*) (1) 

We define a real detectable point source as one which 
yields a positive measurement, so p(P\R) — 1. The prob- 
ability of obtaining a positive measurement depends on 
the false detection rate as well as the real detections. 
Letting F represent a lack of point source (for false de- 
tections), 

p(P) =p(P\R)p(R)+p(P\F)p(F) 

= p(P\R) P (R)+p(P\F)(l-p(R)) (2) 
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Fig. 4. — Best-fit measured distributions of mean surface brightness (left; see text for definition) and point source to host galaxy flux 
ratio (right, expressed in terms of magnitudes — objects with low values of m p — m h have a more dominant point source than those with 
high values). The horizontal bar in each plot shows the range of values encompassing the 10th to 90th percentiles of the distributions. 
Distributions for control galaxies are shown as dashed lines. The m p — distribution for control galaxies comes from fits with a point 
source + galaxy model on inactive galaxies which should not have a real central point source. See text for further discussion. 

probability of having a real detected point source. 

Figure [6] shows the host-to-point source contrast as a 
function of redshift for X-ray and radio AGN. By looking 
at the distribution of points in the y-direction (contrast), 
we see that a substantial fraction (~47%) of X-ray AGN 
fall below our conservative contrast cut, making them 
inconsistent with normal galaxies lacking a point source. 
The radio AGN, however, have a distribution broadly 
consistent with that of normal galaxies, so they do not 
have detectable optical nuclear point sources. Neither 
sample shows a strong trend with redshift. This sug- 
gests that our fitting procedure, whose success (as dis- 
cussed below) depends on the particular distribution of 
light, does not exhibit strong selection effects with red- 
shift over the range considered here. We show later that 
despite the total light output being dominated by stars, 
the nuclear point source can significantly impact mea- 
sured morphologies. 




Fig. 5. — Cumulative distribution of point source to host galaxy 
contrast measured for control galaxies which should lack a real 
point source. Vertical lines show our chosen contrast cuts to sepa- 
rate real point source detections from spurious ones, in the aggres- 
sive (thick dashed, m v — < 3.4) and conservative (thick solid, 
m p — mh < 3.0) cases. We detect point sources brighter than 
the conservative contrast cut in fewer than 7% of control galax- 
ies, so AGN candidate point sources measured to be brighter than 
this contrast level are unlikely to be spurious detections. The thin 
dotted line shows the distribution for galaxies with late spectral 
types only, and the dashed early spectral types; the variation due 
to galaxy type is minimal. 

The unknown probabilities in this equation can be esti- 
mated from our fits and the choice of L. Choosing L = 
3.0, the probability of false detection is p(P\F) ~ 0.07 
based on fits to normal galaxies. The probability of mea- 
suring ni p — mh < 3.0 can be gleaned from our fits to real 
AGN. As shown in Figure H about 30% of AGN fits yield 
a result with m„ — mh < 3.0. Taking p(P) = 0.3, we can 
solve Equation |2] for p(R) : 



P (R) 



p{P)-p{P\F) 0.3-0.07 



p{P\R) - p{P\F) 1-0.07 



: 0.25. 



(3) 



We thus use Equation Q] to find p(R\P) ~ 0.83, so a 
random AGN with m p — < 3.0 has roughly an 83% 



3.3.1. Reliability and Systematics 

We use the simulated AGN images described above to 
constrain systematic uncertainties in our fit results. We 
fit each of the PSF simulations nine times, each time 
using a different model PSF taken from the grids de- 
scribed above. By comparing the resulting distributions 
of best-fit parameters to the known original parameters, 
we can characterize the effect of PSF variations on our 
fit results. Since we intentionally choose incorrect PSFs 
in this test we expect the results to be markedly worse 
than for our fits to real images. Figure [7] shows the dif- 
ferences between input and best-fit parameters for the 
PSF simulations as a function of the difference between 
point source and host galaxy magnitude. As expected, 
our ability to recover accurate parameters with a mis- 
applied PSF is generally better when the host galaxy 
dominates the flux of the entire system. This does not 
hold for the point source magnitude for bulge-dominated 
systems, where the compact galaxy profile can mimic a 
point source. In general, when the point source domi- 
nates the host galaxy light, the uncertainty in the fits 
increases with smaller galaxy radii and more compact 
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Fig. 6. — Point source to host galaxy contrast (expressed in terms of magnitudes) as a function of redshift for X-ray selected AGN 
(left) and radio-selected AGN (right). Error bars are derived from fits to simulated AGN images, as described in the text. The thick solid 
line separates images with real point source detections (m p — m/, < 3.0) from those with likely false detections. While X-ray AGN have 
relatively brighter point sources than radio objects (which are consistent with objects lacking a point source), there is no significant trend 
with redshift. 



profile. The difficulty of recovering parameters for bright 
point sources is likely accentuated by the inaccuracies of 
the model PSF wings, which can become confused with 
light from the host galaxy. Since most of the real AGN 
in our sample have m p — > 0, and we have taken 
care to choose the best PSF, we expect these effects to 
be minimal. 

The recovery simulations help delineate the reliability 
with which we can recover parameters in the presence 
of realistic noise. From the fitting results for these im- 
ages, we determine uncertainty estimates as a function 
of host galaxy mean surface brightness. We fit each im- 
age as we would a real AGN image, using the correct 
PSF. Since the occurrence of galaxies is rare in the ran- 
domly selected noise images, we do not perform mask- 
ing or simultaneous fitting of additional galaxy models. 
We ignore fits that did not converge, and discard fit re- 
sults that yield unphysical parameter values outside the 
boundaries set by the parameter constraints (such cases 
account for ^15% of the fits). The results are shown 
in Figure [8] We clearly recover parameters better for 
brighter host galaxies, with substantial decreases in reli- 
ability of magnitude, radius, and Sersic index for fainter 
objects. Notably, point-source magnitudes are recovered 
to within a tenth of a magnitude in all cases, indicating 
that background noise is a minor problem compared to 
PSF accuracy (characterized in Figure [7]) when measur- 
ing this parameter. 

Our method for assigning realistic uncertainty esti- 
mates to the best-fit parameters of our a ctual AGN 
closely follows that of lHauBler et all (|2007t ). First, we 
calculate the mean surface brightness of the AGN host 



galaxy and take the corresponding standard deviations 
of the parameter distributions shown in Figure [5] For 
any one parameter, we have two uncertainty estimates, 
a n =i and cr n= 4, corresponding to the exponential disk 
and de Vaucouleurs simulations. We linearly interpolate 
between these two values to match the measured value 
of n, using the limiting values instead of extrapolating 
for n > 4 and n < 1. Because our objects are mostly 
host-dominated, we do not attempt to determine how the 
uncertainty estimates vary with the point source bright- 
ness. The additional scatter due to the range in relative 
point source brightness is already folded into the scatter 
in Figure [H] 

The simulation fit results do help us to identify sys- 
tematic effects of a bright nuclear point source. Quali- 
tatively, a bright central PSF causes best-fit values of 
to be systematically low, n to be systematically low, 
to have additional scatter, and m p to have less scatter. 
These effects are pronounced primarily for objects with 
m p — m,h < 0, which excludes the bulk of our sample. 

3.3.2. Morphologies from 2-D Fits 

Table [U gives the results of 394 2-D fits for the AGN 
candidates in our sample, including all objects with ACS 
images in all classes. Initially, 174 fits ran into the fit con- 
straints and we tried again with new parameter guesses 
(see discussion of our fitting methods above). Of the re- 
fits, 74 were successful. After visually assessing the qual- 
ity of the fits, we attempted an additional 74 fits with 
initial parameter guesses adjusted manually, of which 26 
yielded acceptable results. Finally, we checked whether 
2D fits without a point source component were successful 
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Fig. 7. — Results of simulations gauging our ability to recover fit parameters with an imperfect PSF solution. Each panel shows the 
difference between input and recovered parameter vs. input PSF magnitude minus host magnitude. Solid lines show the 10th, 50th 
(median), and 90th percentiles of the recovered minus input parameter difference. The left column shows recovered host magnitude, PSF 
magnitude, radius, and Sersic index for simulated AGN with exponential disks (n = 1), and the right column shows the same for simulations 
with deVaucouleurs profiles (n = 4). The horizontal bar at upper right of each column shows the 10th to 90th percentile range of contrast 
measured for the real AGN in the sample. 



when the fit with a point source failed - eight of those ob- 
jects whose best fit parameters had run into constraints 
with a PSF component did not run into constraints with- 
out. The values in the table are the best best-fit param- 
eters. We assign a flag to each fit, with the following 
values: 0) successful automated fit; 1) fit ran into pa- 
rameter constraints; 2) automated fits ran into param- 
eter constraints, but manual fit did not; 3) fit does not 
include a point source component (because inclusion of 
a point source led to a poor fit); 4) visually assessed as 
a poor fit. We consider those objects with flag values of 
0, 2, and 3 as successful. Table Q] shows the numbers of 
AGN with successful fits in each of our defined classes. 
In the following section we discuss asymmetry and con- 
centration of the AGN, which can be measured even for 
those objects whose 2-D fits have failed. Such objects 
have higher values of asymmetry (by ~40%) but similar 



values of concentration compared to those with success- 
ful fits. It is not surprising that those objects with the 
greatest degree of asymmetry are the most difficult to fit 
with simple galaxy models. 

Figure [5] shows distributions of the best-fit Sersic in- 
dex for X-ray and radio objects and their control sam- 
ples. Class 2 objects have a distribution of Sersic index 
which is statistically indistinguishable from Class 3 ob- 
jects, so we combine them here. For all AGN candidates 
with measured m p — m/, > 3.0 and successful no-PSF 
fits, we replace the Sersic index measured when includ- 
ing a point-source component with that measured with- 
out including a point source. This replacement, which 
accounts for those AGN without a strong nuclear point 
source, ultimately affects the overall distribution of Ser- 
sic index minimally. Control galaxies are separated into 
late and early spectral types using the COSMOS photo- 
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Fig. 8. — Results of simulations gauging our ability to recover fit parameters in the presence of noise. Each panel shows the difference 
between initial and recovered parameters vs. host galaxy mean surface brightness. Solid lines trace the 10th, 50th (median), and 90th 
percentiles of the parameter difference. The left column shows recovered host magnitude, PSF magnitude, host radius, and Sersic index for 
simulated images with exponential disk (n = 1) profiles, and the right column shows the same for simulations with deVaucouleurs profiles 
(n = 4). The horizontal bar in top panel of each column shows the 10th to 90th percentile range of fi measured for the real AGN in the 
sample. 

a fair amount of scrutiny. In particular, our simula- 
tions show (cf. Figure [5]) that a substantial number of 
bulge-dominated AGN will have a recovered Sersic in- 
dex with n < 2.5, the typically used cutoff between disk- 
and bulge-dominated. Quantitatively, <30% of recov- 
ered Sersic indices will incorrectly have n < 2.5. On 
the other hand, 43% (40/93) of the AGN in our mea- 
sured distribution have n < 2.5, indicating a significant 
disk-dominated population. Furthermore, systematically 
low values of measured Sersic index are more prevalent 
for point-source dominated images included in our sim- 
ulations, but most of our X-ray AGN are actually host- 
dominated. Another possible effect emerges from the fact 
that 20% of our X-ray AGN did not yield successful 2- 
D fits. Perhaps the objects with failed fits are exactly 
the ones which would fill in the bulge-dominated por- 
tion of the distribution. Our simulations show, however, 



metric redshift catalog (Mobash er et al.ll2007f ) T phot pa- 
rameter, and we find a good corresponding separation 
of morphologies into disk- and bulge-dominated. Early- 
type control galaxies are clustered around n — 4, al- 
though with significant scatter, and late-type galaxies 
around n = 1. With this division, approximately 60% of 
all the control galaxies have late-type morphologies and 
spectral types, and the remaining 40% have early-type 
morphologies and spectral types. 

X-ray AGN have a Sersic index distribution interme- 
diate between the late- and early- type control galaxies, 
including a broad range of morphologies. When the con- 
trol galaxies are not separated by spectral type, a KS 
test rejects the hypothesis that the X-ray AGN distribu- 
tion is consistent with that of the controls at the 97% 
level. T his result, which conflicts with some previous 
findings (jGrogin et al.ll2005l : iPierce et al.| [2007) . deserves 
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Fig. 9. — Distributions of measured Sersic index for X-ray (left panel) and radio (right panel) AGN, along with matched control galaxies. 
Control galaxies are separated into late (dotted line) and early (dashed line) spectral types based on photometric redshift fits. X-ray AGN 
(shaded, left panel) include Class 2 and 3 objects (with and without narrow lines), and exhibit a morphology distribution intermediate 
between disk- and bulge-dominated. The distribution for radio AGN (shaded, right panel) appears bimodal, with a significant late-type 
morphology component. However, contamination of the radio sample by star-forming and hybrid galaxies, combined with the uncertainties 
in our fits, prevents us from distinguishing the radio AGN morphology distribution from that of early-type galaxies. 

that disk-dominated systems are more likely to fail than 
bulge-dominated ones for all values of the mean surface 
brightness. Finally, we show in the following section that 
measurements of the concentration of these AGN rein- 
forces the trend. 

Like the X-ray AGN, radio AGN have a Sersic in- 
dex distribution spanning a range of morphologies, with 
an apparent bimodality. In this case, 37% (34/92) ob- 
jects have n < 2.5. We emphasiz e, however, that the 
radio AGN classification scheme (|Smolcic et al.l 120081) 
admits —20% contamination from star-forming and hy- 
brid galaxies. This effect might help explain the ap- 
parent bimodality, since the radio objects classified as 
star- forming galaxies have disk-dominated morphologies. 
With this consideration, we cannot rule out that the dis- 
tribution of morphologies for radio AGN is consistent 
with that for early-type control galaxies. For none of the 
samples do we detect any evolution of morphology with 
redshift. 

Because only 18/34 broad- line AGN images were suc- 
cessfully fit (due to difficulties fitting point source- 
dominated objects), we cannot make strong claims re- 
garding their host galaxies' morphologies. The mea- 
sured Sersic index distribution favors disk-dominated 
morphologies, but our simulations indicate that best- 
fit parameters are highly uncertain for such point-source 
dominated objects, especially with an imperfect PSF. We 
used the PSF-only fits to determine which objects have 
resolved host galaxies in the ACS images. To calibrate 
this method, we apply the PSF-subtraction technique on 
-60 stars selected from different COSMOS ACS tiles. 



When we measure the residual flux through an aperture 
after subtraction, we find that residual flux is > 1% of 
the total (pre-subtracted) flux for 16% of the stars. The 
residual flux is > 5% of the total flux for only 3% of the 
stars. Thus, we use a 5% flux cutoff (i.e. if the flux in 
residuals is above 5%, then we claim a host galaxy detec- 
tion, and if the flux in residuals is less than 5%, then we 
do not claim detection) and expect to have false detec- 
tions ~3% of the time. This sort of subtraction technique 
is conservative in the sense that it almost always over- 
subtracts the PSF, yielding a lower limit on the residual 
flux attributed to the host galaxy. Using this 5% tol- 
erance, all of the broad line AGN in our sample have 
resolved host galaxies, though some at a marginal level. 
Table [3] shows our measured limits on PSF magnitude 
and host magnitude for these objects, compared with 
the best- fit quantities taken from Table [2] 

3.4. Asymmetry and Concentration 

We use the asymmetry parameter, A, and concentra- 
tion, C, to further quantify AGN host morphologies. 
These model-independent indices (along with dumpi- 
ness) have been used as "fundamental" properties to 
classify galaxy structure ((Abraham et alj 11994 119961 : 
lConselicell200nh . 

Structures with low spatial frequencies (large scales) 
dominate the asymmetry index, with < 30% o f a galaxy's 
asym metry arising due to star formation (IConselice I 
2003). Therefore large asymmetries serve as good in- 
dicators of recent merger activity, with 50% of nearby 
ULIRGS (expected to be merging systems) showing a 3- 
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sigma deviation fro m the asymmetry trend with colors 
for normal galaxies (|Conselice I l2003h . A conservative 
minimum asymmetry for merging systems is A = 0.35, 
but we do not apply this limit here because we are inter- 
ested only in a difference in asymmetry between active 
and non-active galaxies. 

In this study, we compare asymmetries measured for 
AGN hosts to those measured for control galaxies to de- 
termine whether AGN activity is more li kely to be as- 
sociat ed with mergers and interactions. iGrogin et al.l 
(2005) use similar logic in applying asymmetry measure- 
ments. Because only one filter of ACS data is avail- 
able for most objects in our sample, we pr obe different 
rest w avelengths as a function of redshift. ICapak et al.l 
(|200l . using COSMOS ACS images in both the F814W 
band as well as the F475W band (which was used to im- 
age ~81 square arcminutes), find that asymmetries are 
systematically different when the F475W band samples 
rest frame UV and the F814W band samples rest-frame 
optical light. Measured values are consistent, however, 
when both bands sample optical light or both sample 
UV light. The authors illustrate that the shift in mea- 
sured asymmetry values for the F814W filter occurs near 
z ~1 . where rest- frame U V begins to dominate. Simi- 
larly, [Smcheilit^l] (|2004f ). using Sersic index to classify 
quasar host galaxy morphologies at z ~1, found that 
most objects' optical and UV classifications were the 
same. We therefore expect only small systematic effects 
due to band shifting in the present study. 

We follow th e method of defi ning and measuring asym- 
metry given in lConselice I f|2000f ) . Starting with an image 
cutout with flux distribution S, we rotate the image by 
180 degrees to get a new image, Sisoi an d define asym- 
metry as A = min(X; \S - S 1S0 \ / £ \S\) - A . The sum 
is over all pixels, and we take the minimum asymme- 
try value from a grid of central pixels near the center 
coordinate of the image. Ao is the asymmetry of the 
background, estimated by taking a median of 25 images 
surrounding the primary target. The images used in con- 
structing the background are taken from the same ACS 
tile as the primary, and each has the same size as the pri- 
mary cutout image. For primary targets near the edge 
of a tile, we shift the grid of 25 images so that all images 
fall within the tile's field of view. 

To measure galaxy asymmetry meaningfully for a 
range of redshifts, we must carefully choose the size of 
the image cutout which we rotate and subtract. A simple 
choice would be a constant physical radius, which trans- 
lates directly to an angular size given a chosen cosmol- 
ogy. Since galaxies come in many si zes, perhaps a be tter 
choice is to use a Petrosia n radius dPetrosian 1 1 1976h . or 
a multiple thereof, as in IConselice I (|200oT r The Pet- 
rosian ^-function, rj(r), is defined as the ratio of surface 
brightness at radius r (from the galaxy centroid) to the 
average surface brightness within r. We then denote the 
Petrosian radius as r Vo , the radius at which rj(r) = rjo- 
We choose r/o = 0.2 and measure asymmetries for im- 
age cutouts with this half-width. Because the Petrosian 
radius can give unphysical values for unusual light dis- 
tributions or for images with multiple objects, we set 
a minimum cutout size of 1 arcsecond and a maximum 
cutout size corresponding to a physical radius of 15 hr x 
kpc (~ 3 arcseconds at z = 0.7). These restrictions pre- 
vent unrealistically small (e.g., less than the PSF full 



width at half max) or large choices of cutout radius. 

We measure asymmetries for both the AGN host galax- 
ies and the sample of control galaxies. Because the highly 
symmetric central point source of an AGN biases the 
asymmetry toward low values, we measure asymmetry 
for images with the point source component subtracted. 
We subtract the best-fit model nuclear point source from 
each AGN image. For objects without successful fits, we 
use residual images from our PSF-only-fit subtraction. 
The resulting A distributions for X-ray and radio AGN 
are shown in Figure 1101 including measurements both 
before and after subtraction of the point source. Point 
source subtraction clearly biases the results for the X- 
ray objects toward lower asymmetry, but very little for 
the radio objects. We perform a two-sided K-S test to 
determine whether the AGN and control sample popu- 
lations are consistent with the same underlying distri- 
bution. We find no evidence that AGN have different 
asymmetry distributions from non-active galaxies, with 
K-S test probabilities of 16% and 47% (where a typical 
tolerance of 5% is used to claim the distributions dif- 
fer). The asymmetries for AGN are generally consistent 
with those of non-active galaxies. We find no correlation 
between X-ray luminosity and asymmetry. 

The concentration parameter serves as an alternative 
to the Sersic index to determine whether a galaxy is dom- 
inated by a highly concentrated central bulge component. 
Here we define the concentration as C = 5 log(r> /r< ) , 
with r> = 0.9r no and r< = 0.5r no . Figure fTTI shows the 
relationship between concentration and Sersic index for 
those control galaxies with successful 2D fits. Despite the 
substantial scatter, the overall correlation is clear. Be- 
cause the relationship appears to flatten out for n > 4, 
we fit a parabola to the control galaxy data, with the 
best-fit equation yielding C = 1.262 + 0.244n - 0.0163n 2 . 
For comparison, we plot the X-ray AGN in our sample, 
showing that the trend is comparable. A parabolic fit 
to just the X-ray AGN data yields fit parameters consis- 
tent with those given above. With our definition of C, 
our best-fit parameters show that a delineation between 
late- and early- type of n = 2.5 corresponds to C = 1.8. 

Figure [T2l shows the concentration distributions for X- 
ray and radio AGN. As in our Sersic index analysis, 
we separate control galaxies into early and late spectral 
types. Both X-ray and radio samples include objects 
both with and without emission lines combined into one. 
We also show the distribution of C values before (thin 
solid) and after (shaded) point-source subtraction for the 
X-ray AGN. The presence of the point source signifi- 
cantly biases concentration measurements to high values 
for these objects. X-ray AGN have intermediate values of 
C between those of late and early type control galaxies. 
Radio AGN also include objects with values of C lower 
than that measured for early-type galaxies, but again we 
caution that the radio AGN sample includes substantial 
contamination. These results support those of the Sersic 
index distributions discussed above. 

4. COMPANION GALAXIES 

Kinematically associated neighboring galaxies provide 
evidence for ongoing galaxy interactions. Without de- 
tailed spectral information, we are limited to counting 
neighbors that are within long cylinders seen in projec- 
tion, but with photometric redshift estimates and suffi- 
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Fig. 10. — Asymmetry distributions for Class 2 X-ray and radio AGN samples (shaded) with matched control samples (dashed). Results 
for the AGN before subtraction of the best-fit point source component are shown as a thin solid line. K-S test probabilities that the AGN 
asymmetry populations are drawn from the same distributions as their controls are shown at upper right in each panel. 
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Fig. 11. — Measured concentration, C, and best-fit Sersic index, 
n, for control galaxies (small dots) and X-ray AGN (triangles). 
The solid line shows a best-fit parabola to the control galaxy data 
only. See text for details. 

cient statistics we should be able to discer n signi ficant 
differences among samples. iPatton et al.l (|2000ft pre- 
scribe a detailed method for counting kinematic neigh- 
bors with limited redshift information in flux-limited sur- 
veys. However, our primary interest is not the absolute 
number of companions per galaxy, but rather the frac- 
tion of AGN with a potential physical companion rela- 
tive to that number for normal galaxies. We therefore 
cir cumvent the need fo r the weighting schemes described 
by IPatton et alj (|2000f ) by carefully choosing our normal 
galaxy sample and neighbor criteria. 

We define the maximum projected physical separation 
of close pairs as R™ ax , and use R™ ax = 20, 50, lOO/i" 1 



kpc in three separate trials (see lPatton et aLlfeOOOl 120021) . 
Due to the aperture size of the photometric catalog, 
neighbors within about 1.5 arcseconds (~ 7.5/i _1 kpc 
at z = 0.7) will not be distinguished from the primary 
galaxy, but we expect the asymmetry measure to reflect 
such close companions. We exclude a candidate neigh- 
bor galaxy if its photometric redshift is greater than \a z 
from that of the primary galaxy, where <r z ~ 0.04(1 + z) 
is the uncertainty in th e calibrated photometric redshifts 
(|Mobasher et al.l l2007f) . We must then impose a mini- 
mum lumi nosity for a candid ate neighbor to qualify as 
countable (jPatton et al.ll2000l ). Using the best-fit appar- 
ent magnitude of the AGN host and control galaxies, 
and the COSMOS photometry for the secondary galax- 
ies, we exclude neighbor galaxies more than dm appar- 
ent magnitudes {Iab) fainter than the primary. A choice 
of dm = 2 should restrict our counting to include only 
neighbors which could undergo a relatively major merger 
with the primary. Other choices of dm yield comparable 
results. An alternative to this apparent magnitude limit 
is an absolute magnitude limit tied to M v , the break 
in the V-band luminosity function. We tried counting 
neighbors with this kind of cut as well, with similar re- 
sults. 

Figure [T3] shows the fraction of X-ray and radio AGN 
with a neighboring galaxy within three different search 
radii. Note that the higher search radii are inclusive of 
the lower ones, so the points plotted are not indepen- 
dent. Error bars are estimated for the Poisson case, with 
a n = N 1 / 2 . We do not find significant differences be- 
tween the AGN and their control galaxies' neighbor frac- 
tions. We note that the AGN neighbor fraction appears 
consistently lower than that for control galaxies in the 
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Fig. 12. — Distributions of measured concentration, C, for X-ray and Radio AGN, along with control samples. Control galaxies are 
separated into late (dotted line) and early (dashed line) types. For comparison, we show the concentration distribution of X-ray AGN 
before point-source subtraction (thin solid). Since the radio AGN are all host dominated, we do not show results for measurements before 
point-source subtraction in the right panel. These results mimic those of our 2-D fits, as shown in Figure [9] 
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Fraction of AGN (triangles) with at least one neighbor 
as a function R™ ax . Control 



within a physical separation of R™ ax 
galaxy neighbor fractions are shown as boxes connected by dashed 
lines. 

figure, but changes in the redshift tolerance and mag- 
nitude cut can reverse this effect. No significant trends 
with morphology or luminosity can be discerned with 
the sample size used here. We conclude that AGN are 
no more likely than non-active galaxies to have a near 
neighbor. With sp ectroscopic redsh ifts from the COS- 
MOS VLT survey (jLillv et al.ll2007D . future work will 
more accurately identify kinematic neighbors. 

5. DISCUSSION 

Recent evidence suggests that galaxies hosting AGN 
represent a transitional population, passing from the 
blue cloud to the red sequenc e in galaxy color- 
magnitude space at redshifts z <1 (jJahnke et al.ll2 004b; 
Sanchez et all 120041 ISilverman et all 120081 : iBundv et all 
20071) . AGN host galaxies at both low and high red- 



FlG. 14. — Measured host galaxy Sersic index as a function of X- 
ray luminosity for X-ray AGN. No significant trend is discernible. 

shift are found to be bluer than quiescent elliptic al 
galaxies (iKauffmann et alj l2003t I Jahnke et ail l2004aHbl ; 
ISanchez et al.l l2004h . indicating recent or ongoing star 
formation. Contrary to the findings of some previ- 
ous authors studying morphologies of mod e rate-redshift 
X-ray selected AGN (|Grogin et all 1 120051: iPierce et al l 
2007|) and quasars at low redshifts (jDunlop et al.ll200"3t 



McLeod k, McLeod Il200lh . our results indicate that X- 



ray AGN hosts may be undergoing a morphological tran- 
sition cojicjrr^ent with a transition from blue to red col- 
ors. iPeng et "all (|2006d) came to a similar conclusion 
when studying the host galaxies of gravitationally lensed 
quasars - 30-50% of quasar hosts in their z > 1 sam- 
ple have disk-dominated morphologies . Further, our re- 
sults q ualitatively agree with those of [iCauffmann "et~aT] 
(2003), whose AGN sample had a distribution of con- 
centration index intermediate between that of early- 
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Fig. 15. — Measured Sersic index as a function of host galaxy absolute V magnitude for X-ray and radio AGN. Median error bars for 
the measured quantities are shown in the lower corner of each panel. Magnitude uncertainties are derived solely from uncertainties in the 
2D fits. Radio objects show a slight trend for brighter host galaxies to have more bulge-dominated host galaxies. No significant trend is 
found for the X-ray objects. 

galaxies, more recent studies employ broad-band opti- 
cal colors, optical spectra, X-ray point source luminosi- 
ties, or infrared colors to select AGN. Furthermore, the 
studies mentioned above encompass the full range of 
rcdshift-areal coverage space, from large numbers of low- 
redshift AGN observed in the SDSS, to modest numbers 
of higher-redshift AGN observed in the pencil-beam sur- 
veys like GEMS, to the handful of high-redshift AGN 
hosts observed through strong lenses. COSMOS falls into 
the pencil-beam category, yielding a substantial sample 
size at moderate redshifts. Although the X-ray luminos- 
ity range (42 < log L x < 44) of our X-ray sample spans 
moderate-to-powerful AGN, the bolometric luminosities 
of our radio AGN are difficult to estimate and are likely 
lower. This makes comparisons even within COSMOS a 
challenge. 

We can attempt to examine biases associated with lu- 
minosity here. Figure[T3] shows the measured host galaxy 
Sersic index versus X-ray luminosity for X-ray AGN. We 
find no notable trend, suggesting that earlier studies suf- 
fered little selection bias, although optical quasar selec- 
tion is not directly comparable to X-ray selection. The 
lack of a trend of host morphology with X-ray luminos- 
ity might also have physical implications for the merger 
picture mentioned above. As the gaseous galaxies pass 
by each other then converge, one might naively expect 
the highest X-ray luminosities to coincide with final co- 
alescence into a galactic bulge. Our data do not support 
such a scenario, although we must defer analysis of the 
details to theoretical work. 

We also estimate the host galaxy luminosities based on 
our measurements of the apparent magnitude. To com- 
pute the rest-frame y-band magnitudes, My, we use the 



and late-type galaxies. In the e merging merger pic- 
ture developed throu gh simulations (|Springel et alj|2005t 
iHopkins et al.l [2006D . a morphological transition might 
make sense as we witness evolutionary phases of a 
merger-triggered AGN. As the merger of two gas-rich 
disk galaxies proceeds, accretion onto the central black 
hole would power the AGN as the host galaxi es disrupt 
and fo rm into a bulge-dominated elliptical. iHasingerl 
(2008) provides another possible scenario where bulge- 
dominated galaxies accrete new gas from their cosmic 
surroundings, leading to the build-up of a disk compo- 
nent and subsequent feeding of the previously dormant 
black hole. 

The discrepancies between AGN host morphology 
studies likely arise from a combination of sample selec- 
tion effects and biases in analysis techniques. We have 
already addressed one key bias, that arising due to the 
presence of a nuclear point source in many optical images 
of AGN. For narrow-line AGN, the point source is typi- 
cally much dimmer than the host galaxy, allowing one to 
ignore its effects on total optical flux and color for sta- 
tistical samples. Despite its relative faintness, however, 
the presence of a point source can substantially alter the 
shape of the surface brightness distribution. This leads to 
morphological measurements biased toward higher light 
concentrations and lower asymmetry. 

Selection effects derive from the manifold methods for 
identifying AGN, each of which reveals a different facet 
of AGN activity. The properties of the host galaxies of 
AGN vary with luminosity, black hole mass, and prob- 
ably redshift, so different samples are difficult to com- 
pare. Where early quasar host galaxy studies picked 
out the brightest quasars living in the most massive 
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spectroscopic redshifts and assume spectral energy distri- 
butions for the AGN host galaxies. We choose the rest- 
frame F-band because it shifts into the observed 7-band 
near the median redshift of our sample, and because it 
serves as a convenient reference to the absolute F-band 
magnitudes derived for all galaxies in the COSMOS pho - 
tometric redshift catalog. Following Ho gg et all (|2002j) . 
we calculate the K-corrections by applying filter curves 
for the F8UW filter of HST and the Subaru V filter- 
used for COSMOS observations. We calculate the cor- 
rections with both an ellipti cal galaxy and an Sb galaxy 
template optical SED from lKinnev et all (|1996l ). and we 
display results from the early-type template. Given that 
AGN host galaxies have blue colors and recent star for- 
mation, the true SED lies somewhere between the two 
templates considered here. However, at all redshifts con- 
sidered here, the K-correction differs by <0.2 magnitudes 
between the two templates, so the choice of template 
does not strongly affect the results. Figure [TS] shows 
measured Sersic index versus these derived host galaxy 
absolute magnitudes. We note that for both X-ray AGN 
and radio AGN hosts, the distribution of absolute mag- 
nitudes peaks around My = —22, so these galaxies have 
similar luminosities to M v ~ —2 2 (for z = 1, compu ted 
by starti ng from the local va lue of Brown et al. 2001] and 
following iCapak et al.l[2007l and lSmith et al.ll2005l in al- 
lowing 1 magnitude of passive evolution to z = 1). We 
see a weak trend (with correlation coefficient « —0.2 for 
X-ray and w —0.4 for radio AGN) of morphology with 
host galaxy luminosity in both samples, where brighter 
host galaxies have bulge-dominated morphologies. This 
trend is not redshift-dependent, and may reflect the gen- 
eral galaxy population. 

Although the major merger picture is elegant and en- 
ticing, none of our AGN samples shows enhancement 
of the merger and interaction indic ators applied. This 
roughly agrees wit h previous studies (|Grogin et al.ll2005t 
iPierce et al.ll2007t ). Furthermore, differences in this re- 
sult between subsamples are not significant: radio and 
X-ray AGN candidates all follow the same trends as non- 
active galaxies. These results suggest that major galaxy 
mergers do not play the dominant role in triggering AGN 
activity, with the likely alternatives being minor merg- 
ers and i nteractions, and dynamical instabilities within 
galaxies (|Hasinger Il200l. cf/). 

We caution, however, that the tools we apply here may 
be too blunt to cut to the heart of the question. The key 
uncertainty in drawing conclusions from tests like these 
is the timescale - for any given merger event, how long 
it takes to go from interaction to merger to coalescence 
to relaxation, when and how long an AGN fueling event 
might occur, and how long interaction indicators will be 
observable. While galaxy counts in AGN environments 
may be connected to the likelihood for mergers, they 
serve as an indirect probe at best. By counting neighbors 
we are finding systems that are likely to merge in the 
future, rather than those that have already merged and 
might be in the midst of AGN fueling. 

Like galaxy counts, the morphological measures used 
in this work may not trace galaxy mergers as sensitively 
as necessary to distinguish recently merged systems from 
normal galaxies at moderate redshifts. Certainly we can 
be confident that high-^4 galaxies are undergoing merg- 
ers, but not all recent mergers necessarily have large val- 



ues of asymmetry. As IConselicel (|2003l) suggests, per- 
haps only systems in certain phases of the merger process 
exhibit the large-scale asymmetries to which A is sensi- 
tive. If major mergers are indeed the trigger of AGN, 
the triggering lags the merger in such a way that the 
dominant light distribution of the host galaxy appears 
essentially as relaxed as a normal galaxy. This presents 
difficulties becauses the typical AGN duty cycle of 10 s 
years is much shorter than a typical galaxy's relaxation 
time, and similar to the free-fall time on which violent re- 
laxation is expected to occur. Another possibility is that 
large scale disruptions such as tidal tails are present, but 
at such a low surface brightness that they are too diffi- 
cult to see and measure with such crude techniques at 
higher redshifts. 

Perhaps minor mergers are the answer. These inter- 
actions may disrupt a gaseous galaxy enough to cause 
central inflow of gas onto a black hole without signifi- 
cantly altering the observe d distribution of light. As em- 
phasized in the review bv l.Togeel (f2f)M ). for most AGN 
the amount of fuel available is not the problem per se, 
but rather decreasing the specific angular momentum of 
that fuel by 99.99% to feed the central black hole. Mi- 
nor mergers could potentially disrupt the inner regions 
of a gas-rich galaxy enough to initiate fueling, although 
this is not seen in simulations. Such an interaction would 
be difficu lt to detect using the t echniques described here. 
However. ICanalizo et al.l ()2007f ) find faint shell structures 
indicative of a merger in deep HST images of the host 
galaxy of a low-redshift quasar whose morphology had 
previously been considered quiescent. Detecting similar 
structures in larger samples of AGN could help reveal a 
merger-driven fueling mechanism, but such work will be 
challenging at moderate and high redshifts. 

These considerations make clear the need for more 
work in understanding the relationships among dynam- 
ics, timescales, and observable properties of galaxy merg- 
ers and AGN. As the resolution and scale of simulations 
improves, we expect new constraints on the merger mass 
ratios necessary to trigger black hole accretion to emerge. 
In future work, we hope to develop new techniques for 
detecting and measuring low surface brightness features 
that may betray recent merger activity. 

6. SUMMARY AND CONCLUSION 

We explored the host morphologies and environments 
of AGN in the COSMOS field. Using X-ray- and radio- 
selected AGN candidates with confirmed spectroscopic 
redshifts, we analyzed host galaxy structural properties 
as well as merger indicators, probing the connection be- 
tween AGN activity and galaxy interactions. The follow- 
ing summarizes our main points: 

• The central point source in optical images of X-ray 
selected AGN has substantial impact on measured 
structural parameters such as asymmetry and con- 
centration. Insufficient accounting for the point 
source can lead to systematically low asymmetry 
and systematically high concentration. 

• Full 2-D fits and concentration measurements 
which account for the central point source in X-ray 
AGN show that their host galaxies have a broad 
range of morphologies whose distribution is inter- 
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mediate between the bulge- and disk-dominated 
regimes. 

• Although radio AGN hosts also appear to have a 
wide range of morphologies, contamination by star- 
forming galaxies prevents us from distinguishing 
them from early- type normal galaxies. 

• Measurements of AGN host galaxy asymmetry do 
not differ significantly from those of matched con- 
trol galaxies. 

• Neighbor counts around AGN are indistinguishable 
from those around matched control galaxies using 
photometric redshifts. 

These findings do not support the hypothesis that ma- 
jor mergers drive black hole activity, but they do suggest 
that the host galaxies of AGN at these luminosities may 
be in a state of morphological transition. Future work 
by members of the COSMOS collaboration will address 
this possibility in more detail by examing the colors and 



star formation rates of AGN hosts and their relationships 
with environment (Silverman et al. in preparation). 
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a Positions in Right Ascension and Declination are for optical counterparts to the XMM-Newton X-ray point-source catalog ( Cappclluti ct al. 2007; Brusa ct al. 2007) and the VLA radio source catalog 
iSchinnerer et al .1 [20071 ') . 

b Redshifts are derived from optical spectra from Magellan/IMACS and MMT/Hectospec d Trump et alJl2007ft . 

c Best-fit morphological parameters based on 2-D surface brightness fits arc: host galaxy apparent F814W magnitude (m^), host galaxy effective radius in pixels (r^; 1 pixel — 0.03 arcscconds), host galaxy 
Scrsic index (n), and nuclear point source F814W magnitude (m p ). Asterisks denote best-fit parameters for fits which do not include a nuclear point source componenent. 

Uncertainty estimates for the best-fit parameters (shown in parentheses) are based on fits to our simulated AGN images, and give the 1-cr scatter on each parameter for an AGN with the given mean surface 
brightness. 

c Non-parametric morphological indicators are: Asymmetry (A), Concentration (C), and Petrosian radius {_r VQ ). 

^ fit flag — 0:Successful automated fit, l:fit runs into boundaries, 2:Successful manual fit, 3:Fit better without PSF component, 4:Poor fit (subjective) 

s Spectral classes from [Trump et all (j2007T) . The first letter indicates X-ray ("x") or radio ("r") selected objects, and the numbers indicate optical spectral types — l:Broad line AGN, 2:Narrow lines, 23:Hybrid 
Narrow line/Red galaxy, 3: Red Galaxy 
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Results of PSF-subtraction for broad-line AGN. 
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